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ABSTRACT 


Typical acoustic propagation in shallow water environments is dominated by 
bottom-interacting paths. The effects of rough bottom interfaces and sediment volume 
fluctuations are investigated using model simulations A numerical study of low-frequency 
(- 200 Hz) broadband pulse propagation is presented and several characterizations of 
bottom factors are examined In particular, the vanations of the interface rms roughness 
and a volume fluctuation strength constant on the time resolution of the broadband pulse 
are investigated. It 1s shown that interface roughness is the dominant factor in the 


degradation of pulse resolution. 
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I. INTRODUCTION 


Over the past few years, the interest in littoral, shallow water environments has 
increased significantly. The success of many naval systems in these environments depends 
on our ability to predict and interpret shallow water acoustic propagation The sound 
speed profile in shallow water typically consists of an upper mixed layer below which the 
sound speed decreases to the bottom. Typical shallow-water environments are found on 
the continental shelf with approximate depths of 200 m. The wave propagation 1s 
downward refracting, that is, bottom-interacting paths dominate in long-range propagation 
as depicted in Figures | and Figure 2. Etter (199!) summarized the effect of the sea floor 
as. 


e forward scattermg and reflection loss (but complicated by refraction in the 
bottom); 


e interference and frequency effects; 

¢ attenuation by sediments; 

¢ noise generation at lower frequencies due to seismic activity; 
e backscattering and bottom reverberation 


The purpose of this thesis is to quantify the influence of bottom variability on the 
resolution of pulse transmissions in simple, shallow water environments. 

Rough bottom interfaces can be modeled by assuming some realistic form of the 
spectral shape. The forward scatter of the acoustic field can then be _ treated 
deterministically. Fluctuations of the sediment properties, such as porosity, density, 
permeability, and shear strength, determine the sound scattering from the bottom volume. 
Yamamoto (1995) summarized the relationship between spectral strength and sediment 
type. We used this information to generate random two-dimensional wave-number 
spectra of the velocity variability in order to simulate a realistic ocean sediment volume. 


These techniques are described in more detail in Chapter II. 
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Typical sound speed profile in shallow water 
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Figure 1. Typical sound speed profile in shallow water. 


Ray trace in shallow water 
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Figure 2. Ray trace in shallow water. 
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Wave propagation in an inhomogeneous medium is defined by the wave equation. 
The solution can be obtained by several approximate methods, for example, Wentzel- 
Kramers-Brillouin (WKB) approximation, ray approximation or parabolic equation 
approximation. The last method was used throughout this thesis study. Over the past few 
decades, ocean acoustic models have been developed and improved to analyze measured 
data and simulate sound propagation in real ocean environments. Modeling provides an 
efficient means by which to parametrically investigate the projector (source) radiated 
sound pattern under the influence of environmental conditions. A perfect model is one that 
perfectly represents reality. In practice, such a perfect model would be as complex as the 
problem it 1s attempting to represent. So, there are some assumptions built into every 
computer model This thesis uses a version of the PE model developed at the University of 
Miami known as UMPE (Smith and Tappert, 1994). 

In order to investigate the effects of bottom interface roughness and volume 
fluctuations, we simulate broadband pulse propagation for four different scenarios as 


follows: 
e flat bottom with homogeneous sediment; 
¢ rough bottom interface with homogeneous sediment, 


e flat bottom intertace with realistic sediment sound speed fluctuations, and 


e rough bottom interface with realistic sediment sound speed fluctuations. 


ioe) 





Il. BACKGROUND 


A. THE ACOUSTIC MODEL 


Because propagation ranges are typically much larger than ocean depths, it is 
appropriate to use a cylindrical coordinate system for long range propagation. Customarily 
we define z as depth, r as range outward from the origin, and @ as the azimuthal angle 
(O—>27) defining direction along a radial in the horizontal plane as depicted in Figure 3. In 
general, full 3-D modeling of the propagation is seldom done due to mathematical and 
computational complexity. It is often a good approximation to the real ocean acoustic 
environment to consider it azimuthally invariant. te, to assume that environmental 
parameters such as sound-speed profile, water depth, and bottom composition are 
independent of azimuthal angle to the source. With this assumption, we concentrate the 
analysis in depth versus range. 

As mentioned in the introduction, the PE method is used in this thesis to model 
acoustic propagation. The parabolic equation method was introduced into underwater 
acoustics in the early 1970s by Hardin and Tappert (1973) who devised an efficient 
numerical solution scheme based on fast Fourier transforms. We will derive a PE solution 
in a 2-D constant-density fluid environment where we have assumed azimuthal svmmetry 
and hence no dependence on the ¢ coordinate. 

We begin with the wave equation for a constant-density medium in cylindrical 


coordinates (7,z) independent of ¢ 


| ee 


ae 
va : 
eg 





= 0 (1) 
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We now assume we have a steady state, time-harmonic acoustic field 


plr,2, at) = plr,z)e"™ (2) 


Top View 


Ocean Surface r 





Side View 


Ocean Water Column 


OO ae 


Ocean Bottom 


Figure 3. Cylindrical coordinates in ocean. 


Substituting this into Eq. (1) leads to the Helmholtz equation, 


ep Nae any 


~ —+tkor p=0 , 3 
ae a = Os — (3) 


re, 
where p(fr,z) is the acoustic pressure, 4, =— 1S the reference wave number, 


=) 


( Cc < * » . * 
Wrz) ~ ae) is the index of refraction, c,is the reference sound speed, and c/(r,z) is the 
ay 


acoustic sound speed It ts within c/7,z) that all features of the environment are 
represented 
Tappert (1977) assumed the solution of Eq (3) in terms of an outgoing Hankel 
function with an envelope function ywr.z). 
Daeg ie y(r.c) Hi '(kyr)} (4) 
where y(r,z) 1s assumed to be slowly varying in range. In the far-field (A,r >>1), the 
Hankel function can be approximated by 


z) 
H(kr)= a (5) 
0 





thus an alternative relationship between the acoustic pressure, p(r,z), and the slowly 


modulating envelope function, y(r,z), 1s 


Ry 
plr,z) = Boy wlrzje™ (6) 


This is the standard definition of the so-called "PE field function" y scaled such that at 
oe yv = land P| = P.. Substituting (6) into (3) yields 


a rl es 1 ee 


a Bi ko (ar = Iy = 0 (7) 


ic {} 
a ae ar cs 





Now we introduce the crucial paraxial approximation based on the slowly modulating 


character of y, 
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Oy , ow 
, << 2ik, a (8) 
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We then obtain the tollowing wave equation, 


21k 


~ =e 
a 
——— + 
{} Ae, 


ee i ~ Vy =Ol5 (9) 
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which 1s the standard parabolic equation. This can be rewritten as 


Cy ik, (, 7 
ages 





Or 2x aera 
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(10) 


which is the form first introduced to the underwater acoustics community by Tappert 
(1974) 

Equation (10) is a one-way wave equation which can be solved by a range- 
marching numerical technique. The UMPE model applied in this thesis uses the split-step 
Fourier (SSF) method (Hardin and Tappert, 1973) which solves the parabolic equation by 


Fourier transform techniques. Define the complex Fourier transform pair as 


17d Wad a I (rk Je®= dk. i Tu(r.4.)] ; (11) 


wrk) = x ; w(r.z)e ope I'l y(r,2)| . (12) 


Cee 


where k- is the vertical wave number. We may rewrite Eq. (10) as 


Oy, 
ree +U.,\W (13) 
where 
(re 
ly =a (14) 
DESO 


represents the kinetic energy operator and 


ue == (07 = 1) (15) 


represents the potential energy operator. This representation of the operators as kinetic 
and potential energy operators is especially insightful when one wishes to forin the ray 
equations which have Hamiltonian form (Smith, [992) 


In the PE/SSF method, solutions of y are marched out in range by 


Ar + Ar)=e th yArTop oY op v(r) (16) 
The operator (/,, 1s simply a multiplication operator in z-space but 7,, contains cross 
terms which couple different depths. However, /,, 1s a multiplication operator in vertical 
wave number space. It is therefore desirable to separate the application of each operator, 
One in z-space and one in é-space (here, & 1s the vertical wave number) In other words, 
the PE field function yw 1s specified at some range r initially in the z-domain A 


transformation is made to the A-domain followed by a multiplication of the k-space 


-~ 


ky WT, ryt : ; . : : . 
Operator e " ” where /,,, 1s the corresponding kinetic energy operator in the wave 
number domain The result is then transformed again to the z-domain and is followed hv a 


ee A eae ‘ ~ ; . > . : 
multiplication of the z-space operator e "The final result 1s the field function at 


ei tial 1s, 
w(r + Ar,z) = mS as aie saline F( wr, -))} | (17) 


This is the general algorithm of the PE/SSF method. 

Previously, we have assumed the operators took the forms defined by Eq. (14) and 
(15). These forms are commonly referred to as the “standard PE” and are only one set of a 
number of various operator forms In our simulation, we will use a higher order form 
which was introduced by Thomson and Chapman (1983) referred to as the “wide-angle” 


PE approximation (WAPE) with operators given by 


; l f ie ‘la 

‘ =1—|1-|—} | 18 

Lee) | ‘ | k. } | : (13) 
L S 

Car = -[n(z) = | (19) 


We then need to detine the initial conditions for the PE tield function, y(r= 0,2) It 
is desirable to begin the calculation by specifying the source in the A-domain. Thomson 


and Bohun (1988) have shown that a wide angle source can be defined as 


ea Syaees 
2k \ f ky) 
V0.4) = Ak on o/ ' LK. ) ( Ai) “—" 


where k-ky corresponds to a propagation angle of 90°. 
B. TRAVEL TIME ANALYSIS 


The calculation of the time domain arrival structure at a given range ts an integral 


part of the analysis in this thesis. The pressure field in the time domain is represented by 
plr.z.t)=F'|p.lr.z)|= 2, Paltz) (21) 

where p.(r,z)is a single frequency component. The arrival time structure at some fixed 

range r= can be computed using the complex field p,(R,z) for many frequencies and 


then Fourier transforming to obtain p(X,z,1), the set of complex pressure values in the 


time domain To avoid large transtorm sizes, the UMPE model basebands the solutions 
and shifts the center frequency bin to zero. To reduce the effect of side-lobes, we apply a 


Hanning filter over the bandwidth before transforming to the time domain. Additionally, 


R 


i@ 


“of — e © by simply performing the Fourier analysis on the 


we neglect the phase factor e 
field function y,(R,z). Therefore, the time domain is heterodyned around the value 


1, Re, and arrival times are computed as values of “reduced time” or (1-/,). 
C. ROUGH BOTTOM INTERFACE 


The ocean bottom roughness is assumed to be characterized by a 2-D power law 
spectrum at high wave numbers which can be written as 


{ 


Wl k >> ) =ak? . (22) 





“corr 
where /..,, Is the correlation length of the roughness and / is the spectral exponent. (Note 
that our use of wave number 1s now with respect to range scales of interface roughness. ) 


To give the tull spectrum a realistically smooth structure, we assume the spectral torm 


: HM 
—— (23) 
(1 a eile: | : 
The rms roughness represents the average amplitude of the undulated bottom interface 
variation and the correlation length can be thought of as a scale of self-similarity. The 


normalization factor sis defined in terms of the rms roughness o, i.e., 
20 | W,(k)kdk = 0° (24) 


and therefore 


) 
Ov) 


Ren 


= Ino 7... (25) 


(7) 


a 
Ma 
I 
We only need the 1-D roughness spectrum along the slice in the two dimensional 
profiles (depth z and range r). Along the direction of range propagation in cylindrical 


coordinates, we have 


r K 
ANN te HN OI, 26 


where K is the total horizontal wave number and-k& is the wave number along the 
propagation direction. Substituting Eq (23) and Eq (25) into Eq (26), we obtain 


( B-1} 





1A (9 yen Dope | ca Bear sa) cae (27) 
where 
\ alia 
LOLs 
2( B \f di 2 A a ae 28) 
oe Bi - 
a\2 “(er } ‘ ae 


and / (x) 1s the gamma function. 


To avoid the annoying gamma function, the UMPE model assumes the 1|-D 


spectrum has the form 


W (k= (14 8 pe) | (29) 


corr 


The roughness realization is then simply rescaled by the desired value of o to provide the 


appropriate rms value. 
D. SEDIMENT COMPOSITION 


Yamamoto (1995) summarizes the transformation between three-dimensional 
wave number spectra and_ the velocity variabilities which determine the ocean sediment 
properties This transtormation can be simplified tor our 2-D case. The 2-D compressional 
wave velocity structure in the sediment may be written as 


o(x,z) = ¢,(1+ 6,(x,2)} (30) 


where C(x,z) 1s the velocity at a point (x,z), Cc. is the average velocity at the seafloor and 
6.(x.z) 1s the relative variability. The relative variability is modeled as a zero-mean 


quasistationary random function, 1.e., 
(6,(x,2)) = 0 (31) 
(6.(x,2)6(x'.y"))= [| S(& 4 )xcostk,(x— 0) +k (2 ~ 2) dk dk, (32) 
The 2-D power spectrum of the relative velocity variability 1s denoted by SfA,,4,). By 
definition, the variance (5 2 represents the total energy of the power spectrum 
(=F sk kak ak. (33) 


It is a plausible assumption that the vertical power spectrum may be given by a power-law 


form 
S(k,) = Bk? (34) 
where / is the spectral exponent and B 1s the spectral strength constant. The 2-D power 


spectrum given by Yamamoto 1s 


S(k,.k. 


(B+1), 
= Fark,’ they (35) 


where qs the aspect ratio and #2 must be an integer 


To produce a random realization of sediment volume fluctuations, we assume the 
complex wave number spectrum has a random amplitude and phase. Ae’, for each wave 
number component. The random amplitude of each component of the power spectrum can 
then be obtained by A* =-/n(r) where r is an independent, uniformly distributed random 
variable in the interval (0,1). The random phase @ 1s another independent, uniformly 
distributed random variable in the interval (0,27). To avoid infinite magnitudes at low 
wave numbers, the spectrum is high-pass filtered. The filtered spectrum is then Fourier 
transformed to obtain a realization of the sound speed variability The entire field of 
sediment tluctuations is first computed on a tine vrid and then interpolated to match the 
computational grid of the PE/SSF algonthm The entire scheme can be described as 


follows 


e Step 1. Construct the 2-D wave number field (x, (71), k. (m)) according to 


ica 
ere. (36) 
k (7) =(n— Ak 1 ee ae 
Noe 210 ae 
(Zax — Zoor) (37) 


k.(m)=(m—WAk, = m=1,2,...,.M-M 


bot ? 
where 


ro = maximum propagation range, 


max 
N = FFT size in range direction. 
> = maximum compution depth. 
2. = average bottom depth 
MPF lesize ident direction. 


M, = index of average bottom interface. 


e Step 2. Generate two random numbers series for the random phase and 
amplitude, 


A(n,m) = -In(a,,) . (38) 


d(7,m) ae (39) 
where 


a,b = uniform random number distribution 0 — 1. 


e Step 3. Compute the random phase and amplitude for every mesh point in the 
°2-D realized wave number fteld. S'(A (ak. (111)) le. 


"(Ak (i), &. (i) = S(K, (11).k.(n)) Aime" 1 ao (40) 


e Step 4. Add high-pass filter to normalize the exponent amplitude at low wave 
number 


e Step 5. Compute 2-D FFT to yield relative velocity variability, 
5,(r,2) = 2D FFT] S'(k, (a)k,(m))). 
The velocity structure of the sediment realization is then given by: 


oP. z) =O aC Re(5.(7. z)). (41) 
FE. BEAMFORMING 


In underwater acoustic propagation, much of the information about the 
environment and the relative source/receiver locations 1s contained tn the vertical structure 
of the acoustic field. This is especially true in shallow water environments where the 
vertical structure can be defined in terms of a relatively small number of modes of the 
acoustic waveguide. Because these modes are generally dispersive, the ability to resolve 
the vertical structure of the time arrival at the receiver can lead to estimations of the 
corresponding source location. 

Compared to the directivity of single sensors, a beamformer permits one to listen 
preferentially to wavefronts propagating from one direction over another. In this sense, a 
beamformer implements a spatial filter (see, for example, Defotta ef al/., 1988, Ziomek, 
1995) Consider a line array of equally spaced elements positioned along the y-axis in an 
intinite, homogeneous medium with an undetermined number of remote sources. Let the 


output of an element located at the origin of coordinates due to the /// source be +/(t). 


14 


Under the assumption of plane-wave propagation in the far-field, a source from direction 
4 produces the following sensor outputs 


i 
a (e S| + pee (42) 


G 


where 7 1s the array element index, d 1s the spacing, @, is the angle of arrival from the /th 


source, and c 1s the speed of propagation. Beamforming 1s accomplished by applying 
weights and time delays to the individual element signals, and then coherently summing 
them The beamformer output when viewed in the direction 4, 1s then 


on Ale A= cin i 
Da » vi + et eee (43) 


C 


where the wv, are weights applied to each element’s signal. 
The directional response characteristics of the beamformer can be derived by 


considering a complex, plane-wave cw signal propagating across the array, 
eh) ar (44) 
The source produces the following sensor outputs due to the arrival angle 
e (t) = eaten (45) 


3 


Ze 3, a 
Meie A, = ie — is the wave number of the /// source. The beamformer output is then 
f c 





NY 
h (1) - Ly w nn sin(&, } sin (4, )] o W(6)e" . (46) 


where 
Not 
W(8,)= Siemens tte 7) 
n= 
represents the spatial Fourier transform of the array weights. 
A two-dimensional FFT (spatial and temporal) can be used to implement the 
beamforming operation for an equally spaced line array. For example, consider a line array 


having N-8& equally spaced elements as depicted in Figure 4. The signals from each 


US 


element are sampled simultaneously at rate /, to obtain a multichannel time series x,,,,. M- 
point FFTs are taken along each channel to obtain the temporal frequency response, .Y,,(j). 


where / 18 the frequency bin number, re , 


VG)= diwx,e 7 ON Os ys M =I (48) 


a 


The X,,(/) are represented in the figure as NV horizontal rows of M FFT bins. 
Next, the N cross-channel FFT 1s performed with appropriate weighting along the 
array, 1.€., 


N- oe 
X(j,k) = . w,X,(jJe wv) Re eh | (49) 


n=() 
The result of the N cross-channel FFT series to spatial frequency (or wave number) 
provides the phase difference information in depth for specific temporal frequencies. The 


relationships between the various indices and the variables used in the beamformer output 


are 
f, = ! ae (50) 
M- 
ays = 
"e : S 
Ne (51) 
6 =s) fe 4] (52) 
m= SIN a 


where / and & are the indices in the frequency and wave number domains, respectively. In 
practice, @, 1S determined by interpolating between values of 4; . After the data has been 
interpolated to arrival angle/frequency space, an inverse FFT is taken along each angle bin 


to provide arrival angle as a function of travel time 
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Figure 4. Beamforming example. 


From [DeFotta, 1988] 
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Hi. MODEL SIMULATIONS 


As was described in the introduction, the principal characterisuc of most shallow- 
water propagation is domimated bv bottom-interacting paths Sound transmission in 
shallow water has an optimum trequency of propagation at longer ranges as shown by 
Jensen and Kuperman (1983) Optimum frequency is a general feature of ducted 
propagation in the ocean. It occurs as a result of competing propagation and attenuation 
mechanisms at high and low frequencies. In the high-frequency regime there ts increasing 
volume and scattering loss with increasing frequency. At lower frequencies the situation 1s 
more complicated. With increasing wavelength the efficiency of the waveguide to confine 
sound decreases (the cutoff phenomenon). Hence propagation and attenuation 
mechanisms outside the waveguide (in the seabed) become important. In fact, the 
increased penetration of sound into a lossy seabed with decreasing frequency causes the 
overall attenuation of waterborne sound to increase with decreasing frequency. Thus we 
get high attenuation at both high and low frequencies, while intermediate frequencies have 
the lowest attenuation As Jensen's study shows, the optimum frequency ts in the range 
200-800 Hz for a water depth of 200 m and we will use 200 Hz as our center frequency 

We assume that the ocean bottom sediment ts a fluid-like half-space which means 
that it supports only one type of sound wave, 1.¢., compressional waves. The information 
required for model simulations should include the following depth-dependent material 
properties: the sound speed profile in the water column and compressional wave speed, 
attenuation and the density in the seabed. For all cases, we will assume the water column 
is homogeneous with a constant sound speed. 

This chapter presents the model simulation results of four test cases. Case | 
examines the propagation of a broadband, impulsive source over a flat, homogeneous 
bottom. Case 2 examines the propagation in a similar environment with a rough bottom 
interface. We then examine which parameters, 1.e., rms roughness, correlation length or 


ratio of them, dominates the beam resolution Case 3 is also similar to Case | but with 


realistic bottom volume fluctuations The factors which most strongly influence resolution 
in Case 2 and Case 3 will be combined into Case 4 
The basic range-independent parameters used in the UMPE model are listed in 


Table 1. 















Model Type Wide Angle PE 
Source Function Type Wide Angle Source 
Reference Sound Speed 1500.00 m/s 


Bottom Depth OO) MWe - 
Maximum Computational Depth 400.0 m 
Mesh Points For Positive Depths oye 
Maximum Range hkoOOkm 
Step In Range 2.0m 
Mesh Points In Range SOOO 
Source Depth 100.0 m 
Frequency Bandwidth } 100.0 Hz 
Sound Speed In Sediment 1600 m/s 


Table 1. UMPE Model Parameters 





A. CASE 1: RANGE INDEPENDENT 


We begin with the simple case consisting of a constant speed water column and 
fluid-like homogeneous sediment separated by a_ flat interface The corresponding 
transmission loss field for the center frequency is depicted in Figure 5. We note the 


symmetric constructive interference due to the finite number of trapped, propagating 





ie 0 : 
normal modes By Snell's law, there is a critical angle at arccos =70 4° (relative to 


horizontal) below which no energy is transmitted into the bottom. From the UMPE output 
field with beamforming at 10,000 m, we obtain the angle-time plot shown in Figure 6. As 
expected, this energy arrives primarily near horizontal with an aperture of approximately 
+20°. The steeper-angle energy 220° is rapidly lost by transmission into the seabed. 


Furthermore, the lowest modes arrive first, before the higher modes, as would be expected 


lle, 


Case 1 (Flat bottom) 


Figure 5. Range-depth transmission loss field for flat bottom interface. 
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Figure 6. Arrival 


from an analvsis of the modal dispersion in this environment A standard calculation of the 
number of trapped modes was found to be consistent with the number of arriving modal 
peaks when the range of the transmission loss palette was extended to expose the weakest 


modal arrivals. 


B. CASE 2: ROUGH BOTTOM INTERFACE 


A rough seafloor will cause attenuation of the mean acoustic field propagating in 
the ocean waveguide due to incoherent scattering thereby decreasing the mean signal to 
noise ratio. For individual pulse transmissions. the instantaneous signal becomes smeared 
in time and space due to mode coupling. A rough bottom interface is added to the 
environment of Case | in order to investigate the effect of scattering on the instantaneous 
resolution of the pulse in the water column. 

As described in Chapter II, there are two parameters that control the roughness, 
ie, the rms roughness and the correlation length In order to determine the dominant 
factor. we simulated two sets of comparisons In one, the correlation length 1s fixed at 
1000 m and the rms roughness has values of 2, 5, and 10 m. In the other, the rms 
roughness is fixed at 2 m with correlation lengths of 40, 200, and 400 m. Figure 7-12 
show the beamformed angle-time plots at 10,000 m. We can see that increasing the rms 
roughness and decreasing the correlation length will reduce the resolution but it is 
difficult to see which factor dominates. However, the distortion changes look similar and 
suggest that the ratio of rms roughness and correlation length may be the important 
parameter We simulated several realizations with fixed ratio Sei rms roughness and the 
correlation length to confirm our hypothesis, specifically, (rms roughness) (correlation) = 
|. 500, | 100 and 1:20. The beamformed results are displayed in Figure 13-15 Combined 
with the previous simulations, we find that higher ratios retain higher resolution and lower 
ratios degrade the resolution more severely. Figure 16-18 show single frequency (200 Hz) 
transmission loss fields in depth-range slices for various rough interface realizations. 


Compared to the tlat bottom case, we see that there is more energy penetration into the 
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Figure 7. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 1,000m. 
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Figure 8. Arrival angle-time plot for rough bottom interface with rms roughness = 
Sm and correlation length = 1,000m. 
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Figure 9. Arrival angle-time plot for rough bottom interface with rms roughness = 
10m and correlation length = 1,000m. 
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Figure 10. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 40m. 
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Figure 11. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 200m. 
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Figure 12. Arrival angle-time plot for rough bottom interface with rms roughness = 
2m and correlation length = 400m. 


28 


TL (dB) 





eI = 
=) 
= oO 
nok 
N 
il 
= 2 
= => Ta 
® 5 c© 
ais wn 
S Neue 
cB) 
= = 
i 5 2 
D ay = 
: 7 , 
S = 
“ ae) 
= 0) 
LD AC ° 
Il © 
: N 
N 
A 
: a 
3) 





= = = = = = = 

© © © © © om a) 

2 eh > 8 & 
(sap) o[sue [eALLIY 


Figure 13. Arrival angle-time plot for rough bottom interface with rms roughness = 
Sm and correlation length = 2,500m. 
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Figure 14. Arrival angle-time plot for rough bottom interface with rins roughness = 
2m and correlation length = 200m. 
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Figure 15. Arrival angle-time plot for rough bottom interface with rms roughness = 
Sm and correlation length = 100m. 
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Figure 16. Range-depth transmission loss field for rough bottom interface with rms 


roughness = Sm and correlation length = 2,500m. 
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Figure 17. Range-depth transmission loss field for rough bottom interface with rms 
roughness = 2m and correlation length = 200m. 
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Figure 18. Range-depth transmission loss field for rough bottom interface with rms 
roughness = 5m and correlation length = 100m. 


bottom at longer ranges, and the stable interference structure has broken down into a 


complex pattern. 


C: CASE 3: SEDIMENT VOLUME FLUCTUATIONS 


The scattering due to the velocity variability within the sediment is considered to 
be the major scattering mechanism in many geological fields in shallow water based on 
extensive measurements of acoustic backscattering and sediment cores (Jackson and 
Briggs 1992, Yamamoto 1995). The geology of shallow-water seabeds is complicated 
because of the physical, chemical, and biological activities which took place though the 
long geological time spans. The measurements of the velocity and density variabilities in 
the ocean are beyond this thesis subject. For our simulations, we shall use the sediment 
model described in Chapter If with some typical values extracted from Yamamoto's 
work. 

In this case, we use a flat bottom interface and add the sediment volume 
perturbation. According to Yamamoto's model, the bottom sediment perturbation is 
controlled by three parameters: spectral exponent B, aspect ratio a, and the spectral 
strength constant B. In his conclusions, the spectral exponent f and aspect ratio a are 
fairly constant regardless of the sediment type but the spectral strength constant B varies 
depending on the sediment type and location. We will fix @ and B and vary the values of 
B by choosing three typical spectral strength constants, e.g., 2x10°, 2x107 and 2x10” 
(units of B are m'?). The sediment field for strength constant B = 2x107 m'? is depicted 
in Figure 19. The rms value of the sound speed fluctuations in the sediment for this value 
of B is 22.4 m/s. The beamformed results for the three values of B are shown in Figure 
20-22. The pulse shape does appear to degrade slightly as the spectral strength constant 1s 
increased. Overall, however, the sediment volume perturbation doesn't seem to affect the 


pulse resolution significantly. 
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Figure 19. Sediment field for strength constant B=2x107 m 
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Figure 20. Arrival angle-time plot for sediment volume fluctuation with strength 
constant B=2x107 m’®. 
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Figure 21. Arrival angle-time plot for sediment volume fluctnation with strength 


constant B=2x107 m 
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Figure 22. Arrival angle-time plot for sediment volume fluctuation with strength 


constant B=2x10° m’®. 
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D. CASE 4: ROUGH BOTTOM INTERFACE WITH SEDIMENT VOLUME 
FLUCTUATIONS 


We have seen the effect of rough bottom interfaces and sediment volume 
fluctuations separately. We now need to combine these two factors. From Case 2 and 
Case 3 results, it appears that the effect of a rough bottom interface dominates the 
pulseresolution The strongest sediment volume perturbation strength constant, B = 2x107 
m'" was chosen to produce the largest volume effect. The interface correlation length 
was fixed at 1000 m and the rms roughness was varied from 2 m to 10 m. Comparing 


these results to those cases with no volume fluctuations showed little difference indicating 


that these levels of interface roughness dominate the effects on pulse resolution. 
oe PULSE RESOLUTION ANALYSIS 


We have shown that the effect of the bottom roughness and sediment volume 
fluctuation degrades low frequency pulse resolution in shallow water propagation. In this 
section, we will quantify the effect of rough bottom and sediment volume fluctuation on 
pulse resolution in order to obtain a clearer relationship. The conventional angle resolution 
is defined as an angle width of 1/e peak amplitude. In order to simplify the Case 2 analysis, 
we fixed the correlation length at a value of 1000 m and varied the rms roughness. Figure 
23 shows several slices in reduced times of 1, 1.5, 2, 2.5 and 3 seconds from an arrival 
angle-time plot for rms roughness of 0 (the Case 1 results), 2, 5, 10 and 20 m. Each 
reduced time data point is actually an average of five points spread over a time bin width 
of 0.2 sec, 1.e., the datum for the reduced time of | sec is an average of data at 0.90, 0.95, 
|.00, 1.05, and 1.10 sec. Polynomial fits to the data produce the smooth curves shown in 
the figures We can see a good angle resolution at lower rms roughness values. As the rms 
roughness 1s increased. the resolution 1s found to degrade, particularly for the later arrival 
tunes This ts expected since later arrivals correspond to higher modes which interact more 
strongly with the bottom As an estimation of this effect, we note that the angular pulse 
width appears to double for an rms value of ~10 m corresponding to an rms: correlation 


ratto-of 1 100: 


40 


In Figure 24. slices are taken at arrival angles 0, 5, 10, and 20 deg In this case, 
three points over a width of | deg are averaged to produce each data point shown. There 
iS an anomalous point at an rms of 2 m due to unresolved peaks of two modes. This point 
was not used in the curve fit. The time resolution degrades when the rms roughness 
increases for all arrival angles except 0 deg. This is due to less bottom interactions of the 
horizontal propagation. Furthermore, we notice that the time width does not degrade as 
rapidly as the angular width. None of the rms roughness values appear to double the pulse 
width in time Phe peak minimum transniussion toss also increases with rms roughness as 
shown in Figure 25 by roughly 5 - 10 dB at an rms value of 10 m 

To perform a similar analysis of the Case 3 environment, the spectral exponent and 
aspect ratio were fixed. Obviously, the sediment volume fluctuations show more stable 
results with varying strength constant as depicted in Figure 26-28. The angle and time 
resolution change little with increasing strength constant. The peak (minimum) 
transmission loss also remains stable. 

In Figure 29, the peak minimum transmission loss computed with the combined 
effects of a rough bottom interface and sediment volume fluctuations shows very similar 
results as in Figure 25. It is reasonable to conclude that these levels of bottom roughness 


dominate pulse resolution effects in simple shallow water propagation. 
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Figure 24. Time width versus rms roughness for arrival angle = 0, 5, 10, 15 and 20 


degree. 
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Figure 25. Peak minimum transmission loss versus rms roughness for reduced time 
= 1, 1.5, 2, 2.5 and 3 second. 
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Figure 26. Angle width versus strength constant for reduced time = 1, 1.5, 2, 2.5 and 
3 second. 
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Figure 27. Time width versus strength constant for arrival angle = 0, 5, 10, 1S and 
20 degree. 
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Figure 28. Peak minimum transmission loss versus strength constant for reduced 
time = 1, 1.5, 2, 2.5 and 3 second. 
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Figure 29, Peak minimum transmission loss versus rms roughness for reduced time 
: fi ; Sp 
= 1.1.5, 2, 2.5 and 3 second combine with strength constant = 2x10° m'? , 
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IV. CONCLUSION 


In this thesis. we have presented a numerical investigation of the effects of bottom 
interface roughness and sediment volume fluctuations on the resolution of broadband 
pulse propagation in a simple, shallow water environment. Predictions of the acoustic 
propagation were accomplished with a broadband, wide-angle parabolic equation model. 
Fluctuations of the bottom interface and volume were computed from realistic spectra 
characterizing each structure. The strength, or level of variability, of each was varied 
independently in order to determine the relative influence on pulse resolution 

To be consistent with measurement techniques employed in underwater acoustics 
experiments, a beamforming algorithm was written to discriminate vertical arrival angles 
as a function of arrival time. When no fluctuations were introduced, the solution to the 
range-independent problem was found to agree with a standard normal mode analysis. 
When fluctuations were introduced, the symmetrical arrival time pattern became distorted, 
indicating a degradation in the resolution of the pulse. Such degradations are expected to 
limit the ability of some operational naval systems 

When interface roughness was introduced, the results suggested that the main 
parameter influencing resolution was the ratio of rms roughness to correlation length. In 
other words, a high value of rms roughness does not independently imply a degradation in 
the pulse shape if the correlation length is quite large. Similarly, a relatively small rms 
roughness may still significantly affect pulse resolution if accompanied by a very small 
correlation length In contrast, the largest realistic scales of volume fluctuations had little 
impact on the pulse shape A combination of large volume fluctuations and various scales 
of rms mntertace roughness showed that degradations of pulse resolution are dominated by 
the interface structure. While it may be possible that a different environment (e.g., strongly 
downward refracting water column and upward refracting sediment) would increase the 
influence of sediment heterogeneity, such cases would also tend to increase interactions 


with the bottom interface resulting in similar conclusions. 
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An attempt was made to quantify these results by measuring the l/e - width of the 
pulse at various arrival angles or arrival times as a function of the fluctuation strength. It 
was observed that an rms interface roughness of 10 m with a correlation length of 1000 m 
produced a doubling of the angular width of the arrivals and increased the transmission 
loss by as much as 10 dB. The time width was more stable, however. In none of the cases 
with volume fluctuations was the signal found to degrade appreciably. When the two 
mtluences were combined. the results were nearly identical to the case with no volume 
fluctuations 

Overall, these results suggest that bottom interactions in shallow’ water 
propagation are dominated by interface scatter. The degradation of pulse resolution 
depends on a combination of the rms interface roughness and an associated correlation 
length scale. In order to know a priori the limitations of pulse resolution, it is necessary to 
obtain an accurate statistical representation of the sea floor topography. In contrast, the 
sub-bottom volume fluctuations do not appear to significantly influence the forward 


propagation 
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